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I. INTRODUCTION 

After the surprising result posted recently by the OPERA collaboration Q , Naumov and Naumov have 
posted an interpretation after which " the neutrino advance of time observed in MINOS and OPERA ex- 
periments can be explained in the framework of the standard relativistic quantum theory as a manifestation 
of the large effective transverse size of the eigenmass neutrino wavepacket^^ [2| ■ 

In that interpretation the wave packet of the traveling neutrino is described as the product of a 
transverse and of a longitudinal probability function 0, without any curvature, something which is 
certainly true at the location of the production (i.e., the source at CERN), but something that might be 
questioned after propagation has taken place. 

One consequence of the approximation made in [3| is that one side of the wave packet of a neutrino 
that would be traveling towards the target detector with a small angle, would be propagating with a 
superluminal velocity (see Fig. 1 of Q). 

In this note I try and revisit the approximation made in j2l| and its consequences. 



II. PACKET WAVE PROPAGATION : BACK TO BASICS 



Let's consider the paraxial propagation of a wave packet from CERN to Grand Sasso. For the sake 
of the present note, let's consider the neutrino as a regular particle in the frame of the standard model, 
in particular that it be very respectful of special relativity. The propagation of the wave packet is then 
governed by the Klein-Gordon equation. For a multi-GeV energy, sub-eV/c^ mass particle I will allow 
myself to neglect the mass, and eventually get to Maxwell equation : 
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I will then follow the derivation of paraxial propagation, as can be found in any textbook on lasers 
physics such as Ref. Q. Spin, and therefore polarization, are neglected here. Separating the time and 
space variations Q, £{r,t) — 5(r)e*'^*, we obtain the Helmholtz equation : 
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with uj = kc. In the paraxial approximation, that is the subject of this discussion, solutions are searched 
under the parametrization of the form £{r) — U{p, z) exp (—ikz), where p = \/ + describes the 
polar variable in the transverse plane (Oxy), and U{p,z) is a function that described the transverse 
distribution of the field as a function of z. We have factorized the term exp(— ifcz), which is the Irst 
order phase term of something that propagates along z with wavenumber k. Assuming a slow variation 
of the envelope of the intensity along the axis, we obtain the paraxial wave equation 
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The solutions to this equation that have a Gaussian profile U oc exp {—p^/w^) are 
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and 
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w{z) describes the transverse size of the beam. On axis, the field amphtude varies hke wq/w, that is 
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The intensity is a Lorentzian function of z with HWHM 



Zo = TTWg/A, 
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named the Rayleigh length, or the betatron function at the production location. Far from the source, 
w ^ Xz/ttwo, which allow to define an angular divergence w' = w/z X/t^wq = wq/zq. 

The phase term ^i{kz — + kp^ /2R) contains three contributions : 

• kz describes forward propagation (together with —ujt) along z. 

• —0 describes a phase variation with z that induces a phase jump of tt through the wavepacket waist. 

• the last term kp^/2R with R = z + Zq/z is the term of interest. 
Let's examine them further; 

• Close to the source (z « 0), the third contribution is kp'^z/2zQ, of the order of z/zq, that is extremely 
small. The phase surfaces are therefore asymptotically planes perpendicular to z. 

• Far from the source, (z — )■ oo), neglecting the 2nd, constant, contribution, we get a phase that goes 
like k{z + p^/2R). 

For a phase surface, Lp = tpo, referenced by the phase on axis Lpo = (p{p = 0), the sphere with 
radius R is described close to the axis by (z — Za) = —(1 — cos6)R — —9^R/2 with p = 9R, that is 
(z — Za) — —p'^/2R and (Za = ipo/k). The sphere is therefore described by z + p^/2R ~ Cste. 

Phase surfaces are therefore spheres with radius i?, centered on the source, with R increasing linearly 
with z. 

Far from the source, the wavepacket is therefore described by a paraxial spherical wave with a Gaussian 
transverse distribution. 



III. DISCUSSION 



The wakepacket beam being asymptotically a plane wave right after production on a longitudinal range 
of the order of the Rayleigh length zq, and a spherical wave asymptotically for z ^ zq, let's estimate the 
numerical value of zg. 

• For an Heisenberg-limited wavepacket, zq could be estimated to X/{ttw''^), where the de Broglie 
length A is of the order of 10~^ fm for a 10 GeV beam, and the divergence of the order of (1 km / 
730 km) [2], that is lO'^. Tiny. 

• For a non-Heisenberg-limited wavepacket, zq is obtained simply by the ratio of the size of the source 
(centimeters) to the asymptotic divergence (10~^)[2], that is, tens of meters. 

The waist region close to the source, where the wavepacket beam is asymptotically plane, is very small 
compared to the actual fiight length. 
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IV. CONCLUSION 



Using the usual treatment of the propagation of a wavepacket, as developed for example in the frame of 
laser physics, we obtain an asymptotic description as a spherical wave with a Gaussian transverse profile, 
as soon as it has left the region close to source. 

The quantum mechanics isochronous surface is therefore asymptotically the same as the isochronous 
surface in classical mechanics, and the propagation time simply determined by the distance from the 
production point to the detection point, and the velocity of the traveling thing. 

In Ref. 0], no attempt to compute the curvature of the wavepacket was performed, and therefore no 
curvature was found, which lead unavoidably to the claim of a superluminal effect. 
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